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Kinetic equations are introduced for the transition-metal nanocluster nucleation and growth mech- 
anism, as proposed by Watzky and Finke. Equations of this type take the form of Smoluchowski 
coagulation equations supplemented with the terms responsible for the chemical reactions. In the 
absence of coagulation, we find complete analytical solutions of the model equations for the auto- 
catalytic rate constant both proportional to the cluster mass, and the mass-independent one. In the 
former case, £fc = s&(£i) oc £i/k was obtained, while in the latter, the functional form of Sfc(£i) is 
more complicated. In both cases, = /i M (Af M (t)) is a function of the moments of the mass dis- 

tribution. Both functions, Sk(^i) and ft M (M M ), depend on the assumed mechanism of autocatalytic 
growth and monomer production, and not on other chemical reactions present in a system. 
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Nucleation and growth phenomena, resulting in occur- 
rence of a new phase from a homogeneous host phase, 
are ubiquitous in nature [HQ- In many cases, apart 
from processes of coagulation and fragmentation, chem- 
ical reactions are also present in a system. Frequently, 
the chemical reactions account for phenomena studied by 
polymer and colloidal science. 

Nucleation and subsequent growth of metal nanoclus- 
ters in aqueous solution is a topic of considerable cur- 
rent interest, since solution route synthesis is one of the 
most convenient methods of producing transition-metal 
nanoparticles However, applicability of this method 
depends on the ability to control size and shape of the 
produced nanoparticles, which determine their unique 
optical, electronic and catalytic properties. For that rea- 
son, theoretical models capable of predicting the cluster 
size distribution, as well as its dependence on the ex- 
perimentally controllable parameters of the system, are 
required. 

A mechanism of transition-metal colloidal nanoparti- 
cle formation has been proposed by Watzky and Finke 
[3], cf. The WF mechanism consists of (i) slow 

monomer, i.e., the zerovalent transition-metal atom (Bi) 
production due to reduction reaction A — > Bi of a metal 
precursor (A) , usually a transition- metal complex coordi- 
nation compound, (ii) fast autocatalytic reduction reac- 
tion A+Bi — » Bj_|_i taking place on the surface of growing 
metal nanoparticles, consisting of i atoms (Bj), and (iii) 
process of coagulation Bj + Bj ^i±j > reversible or oth- 
erwise. In the original WF scheme [4|, step (iii) had not 
been considered. Irreversible coagulation was first intro- 
duced in 

For the transition metals in which higher oxidation 
states are present, e.g. Au, at least one additional pre- 
liminary step of the form (iv) P — > A is needed [fj-ll3| ■ 

As an excess of the reducing agent is usually used in 
the reactions (i), (ii), (iv), its concentration is fairly time- 
independent, and all chemical reactions may be treated 



as irreversible. Furthermore, (i) and (iv) may be treated 
as reactions of pseudo-first order, while (ii) as a reaction 
of pseudo-second order. 

The WF mechanism is applicable to other systems. 
Certain cases of transition metal oxides or sulfides nan- 
ocluster formation, or polymerization phenomena of var- 
ious kind, including protein aggregation Q, are well de- 
scribed by an effective model defined by (i)-(iv), even if 
the actual mechanism of nucleation and growth in these 
systems is more complicated. 

However, the kinetic equations corresponding to step 
(iii), proposed in 0-01' cannot be regarded as fully sat- 
isfactory. In particular, the original approach does not 
allow to predict the cluster size distribution. Thus, in 
this Rapid Communication proper kinetic rate equations 
for the WF mechanism are introduced. 

Kinetic rate equations. In order to describe kinetics of 
step (iii), the approach proposed by Smoluchowski (hH . 
15| . a standard and widely-used tool for description of 
various coagulation phenomena, is reintroduced (l6l - [29| . 
In the present case, concentrations of P (c„), A (c Q ), and 
Bj, (£i), i G N, are the state variables. Smoluchowski 
coagulation equations have to be supplemented with the 
terms related to steps (i), (ii), and (iv). In effect, the 
time evolution of the system is given by the following 
rate equations 
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where k^ = Cpfc^, fc a = c p fc a , and i?i = c p Ri are observ- 
able reaction rate constants for steps (iv), (i), and (ii), 
respectively, while k n , k a , and Ri are the correspond- 
ing bare rate constants, and c p is a concentration of the 
reducing agent, assumed here to be time-independent. 
Kij = Kyi and Fy = Fji denote coagulation and frag- 
mentation kernels. The first sum in (|4|) is restricted to 
i + j = k. 

Initial conditions of interest here are c n (0) = cq, 
c Q (0) = d , with c + d = <?o 0, and &(0) = 0, 
i 6 N. The solution for c n (t) is immediate: c^(t) = 
Co exp(— k v t). The original WF model corresponds to 
c - 0. 

It follows that c v (t) + c a (t) + J2"jLijt.j(t) — (mass 
conservation constraint). Therefore, the quantity 



(5) 



is conserved during the time evolution. 

Apart from the assumption of constant concentration 
of the reducing agent, c p (t) = c p (0), (.41), it was also 
assumed that no source term (no injection mechanism) 
for either P or A species is present (A2); that both the 
autocatalytic P to A reduction reaction, P + B.; — > A + 
Bj (A3), as well as the disproportionation reaction P + 
Bi < — > 2A (.44) may be neglected. It was also assumed 
that neither K^, nor Fij depend on concentration of c p , 
or c a (45), and finally, that chemical species P, A, 
and reducing agent (R) do not form clusters (46). Any 
of the above assumptions may be abandoned, leading to 
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a generalization of the model defined by Eqs. ©-((4]) 

If R\ = 0, for all i £ N, D and (@| become the stan- 
dard Smoluchowski equations, with the monomer source 
term k a c a (t) for k a 7^ 0. However, for Ri such re- 
duction is no longer possible, and Eqs. (JTJ)— with their 
generalizations arc members of a wider class of 'reaction- 
aggregation' equations. Few models of this type have 
been found in the literature [26l-28|. 

Method of moments. To analyze properties of Eqs. 
©-(jj), it is convenient to apply the standard method of 
moments. The fi-th moment of the cluster mass distri- 
bution is defined as M p (t) = Y^jLi^^ji^)- From Eqs. 
([3]) and ([4]) we obtain 



j=l p,q p=2 



where S$ = S$ = §((p + ?)" - p" - q»)K pq , T P M 
Eti (*" ~ lP»)Fi,P-i and = [(j + 1)" - j»]R r 
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From Eq. (|6]) we see, first, that all moments grow due 
to the monomer production (i), note, k a c a > 0. Sec- 
ondly, for fx, = 0, = 0, because the total cluster 
concentration M is not affected by autocatalytic reac- 



tion (ii). For fx = 1, S t 



(i) 



(i) 



0, i.e., the total 



cluster mass is not changed by coagulation or fragmen- 
tation. Thus, for fx = 1, stationary solution of © exists 
for k a > 0. Indeed, for (k a + J2j=i > 0, from 

Ml = follows c a = 0. Since lim^oo c v (t) = 0, Eq. flSJ) 
gives lim^oo M\(t) = Mi = qo, as expected for an irre- 
versible reaction. Additionally, for Kij = Fij = 0, when 

S P q^ = T~}^ = 0, similar procedure proves that all M^ 
approach stationary value. 

In order to obtain tractable system of time-evolution 
equations for M p (t) and c Q , a restriction is imposed on 
fi and so that fi € NU0, Fij = 0, and Ri and Kij 
given by 

Kij = k + + j) + K 2 ij, Ri = d R i + b R , (7) 

where a R , b R , k , K\, k 2 are arbitrary non-negative co- 
efficients. Under these simplifying assumptions, Eqs. © 
and © for fi = 0, 1, 2 assume the form 



M 

Mi 
M 2 



1 



a R Mic a - b R M c a , (8) 



KqMq - KiMqMi 



(9) 
(10) 



+ [2a R M 2 + (a R + 2b R )M 1 + b R M ] c a 
K Ml + 2k 1 M 1 M 2 + k 2 M|. (11) 



The initial conditions, Vk : ^ fe (0) = give M p (0) = 0, 
V/x, while variables c a and Mi are not independent, since 
they obey the constraint (HJ) for Mi = XwliJ'CjW- 

Explicit form of the corresponding evolution equations 
for higher moments (fx > 2) may be easily found using 
Eq. (|6|). However, higher moments are not needed to 
calculate two basic characteristics of the cluster mass dis- 
tribution, i.e., mean cluster size, (£), and variance, cr 2 (i), 
as given by 

Mi o / \ M 2 / Mi \ 2 , . 

j4&sence o/ coagulation. The case of negligible coagula- 
tion and fragmentation is at the center of our interest. 
In colloidal systems, experimentally, this is achieved by 
addition of a stabilizing agent, e.g. the polyvinyl alcohol 
(PVA) or polyvinylpyrrolidone (PVP), which inhibits 
coagulation, affecting chemical reactions to a lesser ex- 
tent [U [l3[. Models of autocatalytic reaction without 
coagulation may also provide an adequate description for 
other systems, e.g. simple chain polymers. 

For Kij = F tj = 0, Eq. g]) divided by Eq. ([3j gives 



d£k ■Rfc-l&B-l — Rk£,k 

d£i k a - Aid 
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Eqs. (fT3"|) for k = 2, 3 . . . form a system of coupled linear 
differential equations, allowing to determine each as a 
function of £i, = Sfc(£i). If coagulation and fragmenta- 
tion processes could be neglected, Eqs. (fT3"j) . therefore all 
Sfc(£i) functions, for given initial conditions, depend only 
on fc a and Ri parameters, and not on chemical reactions 
and physical process, e.g. injection of P and A substrates, 
unless they involve Bi clusters. In particular, Sfc(£i), 
hence the stationary values of = \im t ^oo £k(t) do 

not depend on assumptions .A1-.A3. 

The above holds true for any equation derived by di- 
viding Eq. or one of Eqs. ©-(HI]), by Eq. © or 
(fTO]) . which allow expressing £i and M^ as the functions 
of M or Mi only, sec Eqs. ((15)1. ([26]). and ((34]) below. 

As a consequence, all special cases or generalizations of 
the model defined by Eqs. CO)-© with the monomer pro- 
duction provided solely by (i), and with identical mech- 
anism of autocatalytic reaction (ii) belong to the same 
universality class. 

Two special cases of the reaction kernel (0 are ana- 
lyzed below: the linear kernel, proportional to the clus- 
ter mass, Ri oc i, b R = 0; and the size-independent one, 
a R = 0. 

Linear reaction kernel. For kq = K\ = n 2 = 0, b R = 0, 
and clr ^ 0, Eqs. (J]), (©, and (fit)]) become identical to 



the rate equations analyzed in Refs. (12|,ll3j (co ^ 0) and 
HHS H| (co = 0). These equations have the same form 
for arbitrary choice of Kij and Fij, hence the reaction 
rate of (ii) , proportional to the total mass of the clusters 
(Mi), does not depend on presence of the coagulation or 
fragmentation. 

c a may be eliminated in favor of c v (t) and Mi = x 
using Eq. §5$, which yields the form of Eqs. (fT 



x = (k a + a R x) (f(t) - x) , 



(14) 



where f(t) = q — c Tr (i) = d + c (l — exp(— k^t)) and 
< x < qo- The following substitution: u — k a + d R x 
transforms (] 14[) into Bernoulli-type equation, which gives 



x(t) 



1 

an 



1 



t 
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where $(t) = Ci(exp(-M) - 1) + C 2 t, C 2 = d R q + 
k a , and C\ — codn/k^ = c§k±jk\. For the original WF 
model, where c = 0, C\ = 0, and f(t) = d ^ 0, lfl"5)) 
reduces to 



Xafi(t) = - [do + ^ ■ J^ e («ndo + k a )t + l \ 

y a-R J \ d a R J 

(16) 

where x a p(fi) = and lim^oo x a p {t) = d 0,E3. How- 
ever, for co ^ 0, x(i) cannot be expressed as a combina- 
tion of finite number of elementary functions. Still Eq. 
(|15p may be given a more convenient form x(t) — x(z(t)), 



where z(t) = C\ exp(— k^t) and xi z ) is given by 



X(z) = H 
a R 



1 F 1 (l- 7 ;2- 7 ;C 1 ) 
^(1-7) 



,7-1 



7-1 



iFi(l-7;2- 7 ;z) \ _ fc„ 
fc„(l- 7 ) J a R 



(17) 



Above, 7=1 + C 2 /k- K and iFi(a;6;z) denote confluent 
hypergeometric function |30j . 

Since in this case, there is an analytic solution ((17)) for 
Mi(t) = x(t), it is convenient to express both M and 
M2 as the functions of Mi. 

Eq. (O divided by Eq. (10]) gives 



cZM 



1 



(18) 



(19) 



(20) 



Using Eqs. ((19)) and ((20)) . explicit formulas for mean 
cluster size and variance (|12p may be easily derived. 
Solving Eqs. (fl3)) for i?i = a^i gives 



1 Qo ( u 



dMi k a + ~a R M 1 1+^Mi' 
where w = qod R /k a = qoa R /k a . Eq. (fl8| gives 



M (Mi) = — In ( 1 + —Mi 
\ qo 



In parallel, dividing Eq. ((TT]) by Eq. ((TU)) gives 



M 2 (Mi) = Mi 1 + — Mi 
9o 



&(*) = * (o) (&(*)) = — > k>l. (21) 

« W \q J 

In order to find £1 (t), M (fT9)) is equated with X^ili £j (*)■ 
Employing ((2"T)) and the identity 



fe=i 



= In 1 + 

x \ 1 — x 



(22) 



where x = Q ■ (O + 1) 1 and O = ujMi/qo, gives 



'/o 



Mi(t)" 



Combining Eqs. ((21)) and (|23|) . the following is found 

fe 



t^OD fc UJ \ LU + 1 



Finally, for c = 0, Eqs. (JTBJ), (JSTJ and ((23]) yield 



4fc k uj \uj + 1 



I _ e -k a (l+v)t 



(24) 



(25) 



Size-independent reaction kernel. For colloidal sys- 
tems, Rj given by Eq. (J7J), for 5^ = and 67^ 7^ 0, 
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provides a lower bound for any realistic functional form 
of Rj , which is expected to be a non-decreasing function 
of j. Moreover, in the absence of coagulation, i.e., for 
K o = K i = K 2 = 0, Eqs. (|5|)- ([TT|) may be regarded as a 
simple model of linear polymer growth [32| , cf. pjj . In 
this case, to solve Eqs. ©-(HI]), Eq. HTDJ) is divided by 
Eq. ©, giving 



obtain explicit form of £1 (t) , Eq. (jT3|) is divided by Eq. 
©. This yields 



Mo 



dM 



5o 



(34) 



The solution of Eq. (|34j) is given by 



dM x 



= 1 



b R 

r Mo = 1 



-Mo, 



dM k a qo 

where uj = qobu/k a = qobu/k a . Eq. (|26p gives 

Mi (Mo) = ^M 2 + M . 
250 



(26) 



(27) 



Applying a parallel procedure to Eqs. (fTTj) and ©, with 
the use of Mi (Mo) (|2"T|) . and solving the equation gives 



M 2 (M ) 



Mo + ^M 2 
^5o 



3al 



(28) 



Using Eqs. ([27]) . (fit)]) and the constraint ([5]) in order to 
eliminate c a , yields the following equation 



M 



k a [qo- Cn{t) - Mo 
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2^ J 
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For Co 7^ 0, it was impossible to find an analytical 
solution of Eq. (|29|) . However, the stationary solution of 
Eq. (|2U]) may easily be found. Combining Eq. (|27l) and 
the limit lim4_ ) . 00 Mi = qo, gives 



lim 



M„(t) + ^M 2 (i) 
2<7o 



= Mo + — M 2 = go. (30) 

2^0 



Eq. ([3D]) implies 



Eq. 



M = ^(vTT2^-l). (31) 

UJ 

may also be easily solved for Co = 0, giving 
(r? - 1) fl - e-^'>' N 



M„(t) = ^ 

UJ 



(32) 



where 77 = VI + 2w. Mo (i) (]32[) has the following proper- 
ties: M o (0) = 0, M (t)_< M , where M is given by ([31]) 
and lim t _ >00 Ma(i) = Mo, as might have been expected. 
Eqs. (|27]) . (|28]t. and (J3J} may be used to calculate the 
time-dependence of (i) p and cr 2 (i) (fT2]l . 
For Ri = bft, solutions of Eqs. (fT3|) read 



5(i 



ft = s w (6) = -U + "(fi)Wfc-i 

UJ 



In 



(33) 

where = ~£i — 1 and Wi(x) = ^}=o x ^ Ifii which 

may be easily verified (note, Wk-i(x) = W' k {x)). To 



M (t)) = £(l- 



(35) 



where for co = 0, Mo(t) is given by ([32]) . 

Using Eqs. ([3T]) , ([33]) . and ([35]) . stationary value of 
each £fc may easily be found 



£ fe = lim Cfc = — 

t— too uj 



l- e -(v-i) Wk _ l{j] _ 1 ) 



(36) 



Finally, for the cluster size-independent reaction ker- 
nel, coagulation process, when present, decrease the rate 
of the autocatalytic reaction (ii). This becomes intu- 
itively clear, when the present model is applied to de- 
scribe growth of linear polymers with two active reaction 
sites at the ends of the polymer chain, since each coag- 
ulation event reduces the number of reaction sites by a 
factor two. Clearly, the influence of fragmentation pro- 
cesses is exactly opposite. 

Summary and discussion. In this paper, time- 
evolution rate equations for the model of transition- 
metal nanoclustcr formation, as proposed by Watzky and 
Finkc, have been introduced. The equations introduced 
constitute a natural generalization of both the Smolu- 
chowski coagulation equations, and the rate equations, 
describing the kinetics of monomer production, autocat- 
alytic nanoparticle surface growth, and other chemical 
reactions present in the system. 

In the absence of coagulation and fragmentation, exact 
solutions of the model equations have been found for the 
autocatalytic rate constant (reaction kernel) proportional 
to cluster mass, Ri oc i, as well as for the cluster-size 
independent one, Ri = const. 

Secondly, it was demonstrated that the functional de- 
pendence of the fc-atom cluster concentrations, on £1, 
given by £ k = s(6-)i an d £k = limt-yoo &(*)> arc com- 
pletely determined by the assumed model mechanism of 
monomer production and autocatalytic growth, but do 
not depend on any other chemical reaction. However, 
this does not hold for the £fc(t) functions. 

In conclusion, kinctically inequivalent generalizations 
of the present model are divided into universality classes, 
whereby two such models belong to the same class when 
they yield the same s(£i) function. 
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